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ABSTRACT With increasing global concerns toward environmental protection and_ sustainable 
development, green supply chain management (GSCM) has drawn much attention from academicians and 
practitioners. Selecting an optimal green supplier is a critical part of GSCM, which can be viewed as a kind 
of multi-criteria decision making (MCDM) problem. To derive the best result, large group of decision 
makers are often involved in the green supplier selection nowadays. Besides, decision makers tend to 
express their evaluations utilizing uncertain linguistic terms due to the vagueness of human thinking. Hence, 
this paper aims to propose a hybrid MCDM approach for green supplier selection within the large group 
setting. More concretely, interval-valued intuitionistic uncertain linguistic sets (IVIULSs) are applied for 
assessing the performance of green suppliers concerning each criterion. Ant colony algorithm is utilized to 
cluster decision makers into several subgroups. The linear programming technique for multidimensional 
analysis of preference (LINMAP) is adopted for the determination of the optimal weights of criteria 
objectively. Finally, an extended MULTIMOORA approach is utilized to generate the ranking of alternative 
suppliers. The practicality and usefulness of the developed large group green supplier selection framework 
is illustrated using an empirical example of a real estate company. 


INDEX TERMS Supply chain management, Green supplier selection, Ant colony algorithm, Interval- 
valued intuitionistic uncertain linguistic set, MULTIMOORA method, LINMAP approach. 


Il. INTRODUCTION 

Because of growing public concerns over the environmental 
issues and government regulations toward sustainable 
development, green supply chain management (GSCM) has 
become popular in production operation management of 
modern enterprises [1]. By integrating environmental 
concerns into supply chain practices, GSCM is considered 
as a promising approach to improve the commercial benefit 
and environmental performance of organizations 
simultaneously. The major aim of GSCM is to decrease 
environmental pollution and eliminate waste during the 
process of purchasing, manufacturing, distributing and 
selling products [2]. As a result, there are a lot of programs 
for companies to implement GSCM, including green design, 
green production, green transportation and green marketing 
[3]. But an organization’s environmental performance 
depends on not only its own green efforts but also its 
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providers’ green practices. Thus, selecting the optimal 
green supplier is a vital part of the GSCM for firms to 
develop sustainability [4-6]. 

Green supplier selection determines the best supplier 
which is capable of providing the buyer with high quality, 
low cost, quick return and good environmental performance 
simultaneously [7]. It plays a vital role in maintaining the 
competitive advantages of a company [8]. An appropriate 
green supplier makes a great difference in enhancing 
quality of end products and satisfaction degree of customers. 
The application of green supplier selection can be presented 
in the situation of multiple suppliers throughout a product’s 
life-cycle [9]. Consequently, it is significant for 
organizations to choose the environmentally, socially and 
economically powerful suppliers. Moreover, a large group 
of decision makers from different backgrounds should be 
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involved in the green supplier selection process due to the 
increasing complexity of the problem. 

In previous studies, many scholars have applied fuzzy set 
theory [10] and intuitionistic fuzzy sets (IFSs) [11] to cope 
with the vague evaluations of decision makers and improve 
the effectiveness of green supplier selection. However, 
when transforming linguistic assessments of alternatives 
into fuzzy numbers, the information of decision makers’ 
subjective judgements may be lost and distorted [12]. As 
for IFSs, the membership and non-membership degrees are 
signified by exact values, which cannot well handle the 
uncertainty and fuzziness of assessment information. 
Recently, Liu [13] proposed the concept of interval-valued 
intuitionistic uncertain linguistic sets (IVIULSs), which are 
the combination of uncertain linguistic variables and 
interval-valued intuitionistic fuzzy sets (IVIFSs). The 
linguistic variable, membership degree and nonmembership 
degree of each element in an IVIULS are denoted by 
interval values rather than crisp values, which can better 
capture the ambiguity and uncertainty of evaluation 
information. Owing to its characteristics and merits, the 
IVIULS theory has been utilized in different fields [14-16]. 
Therefore, it is natural to use the IVIULSs to evaluate the 


green performance of alternative suppliers on each criterion. 


On the other hand, selecting the best-fit green supplier is 
often viewed as a multi-criteria decision making (MCDM) 
problem and a variety of MCDM approaches have been 
used for green supplier selection in recent years [1, 3, 4, 6]. 
The multi-objective optimization by a ratio analysis plus 
full multiplicative form (MULTIMOORA) method is a 
distinctive MCDM technique proposed by Brauers and 
Zavadskas [17]. It includes three parts as follows: the ratio 
system, the reference point and the full mutiplicative form 
[18]. The final ranking of alternatives is made in 
accordance with the dominance theory [19]. Comparing 
with other MCDM _ methods, the merits of the 
MULTIMOORA approach are that: (1) It can make the 
ranking result more accurate by aggregating the three basic 
parts; (2) It can effectively solve complicated MCDM 
problems with numerous alternatives and criteria; (3) Its 
calculation process is easily comprehensible and the result 
can be obtained rapidly. Hence, it is desirable to apply the 
MULTIMOORA method to address green supplier 
selection problems. 

Against the above discussions, we propose a hybrid 
MCDM model based on IVIULSs and MULTIMOORA 
method in this paper for evaluating and ranking green 
suppliers under large group environment. We utilize ant 
colony algorithm to cluster decision makers, an extended 
LINMAP approach to calculate the objective weights of 
criteria, and an improved MULTIMOORA technique to 
rank alternative green suppliers. The presented large group 
green supplier selection approach is able to reflect the 
vagueness and ambiguity of decision makers’ judgements 
and acquire an accurate ranking result of green suppliers. 
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The remaining sections of this paper are arranged as 
follows: In Section II, we review the existing green supplier 
selection approaches and the applications of the 
MULTIMOORA technique. The basic concepts and 
definitions regarding to IVIULSs are introduced in Section 
II. In Section IV, we put forward a hybrid MCDM model 
using IVIULSs and MULTIMOORA method for large 
group green supplier selection. In Section V, a case of a real 
estate company in China is given to demonstrate the 
proposed approach. The last section presents concluding 
remarks and future research recommendations. 


Il. LITERATURE REVIEW 


A. GREEN SUPPLIER SELECTION METHODS 

In recent decade, a lot of green supplier selection methods 
based on MCDM have been developed, which can be 
classified into individual methods and hybrid methods [20]. 
As for individual methods, Sanayei et al. [21] utilized fuzzy 
VIKOR (VIsekriterijumska optimizacija i KOmpromisno 
Resenje) method to solve the supplier selection problem. 
Chen and Zou [22] applied intuitionistic fuzzy grey relational 
analysis (GRA) to find the best-fit supplier from the 
perspective of risk aversion. Cebi and Otay [23] used the best 
worst method (BWM) to choose the optimum green supplier 
for an edible oil company. You et al. [24] used a modified 
VIKOR method for green supplier selection with interval 2- 
tuple linguistic information. Keshavarz Ghorabaee et al. [25] 
proposed an improved complex proportional assessment 
(COPRAS) technique to select the optimum supplier in the 
interval type-2 fuzzy context. 

In addition, many hybrid MCDM models have been 
utilized to deal with green supplier selection problems. For 
instance, Gary Graham et al. [26] offered a new model by 
integrating analytic hierarchy process (AHP) and TOPSIS for 
selecting the optimal green supplier. Wang et al. [27] 
proposed a hybrid MCDM model based on cloud model and 
QUALIFLEX method to evaluate the green performance of 
suppliers. In [28], intuitionistic fuzzy AHP and intuitionistic 
fuzzy axiomatic design were integrated to handle the green 
supplier selection problem for an international sporting goods 
firm. Tsui [29] developed an innovative MCDM model by 
combining DEMATEL-based analytic network process 
(DANP) and PROMETHEE to select the best green supplier 
in the thin film transistor liquid crystal display industry. Qin 
et al. [1] presented the TODIM technique based on prospect 
theory for the selection of green suppliers under interval 
type-2 fuzzy setting. Yazdani [30] designed a green supplier 
selection model, where DEMATEL method was utilized to 
handle the inter-relationships between customer requirements, 
quality function development (QFD) technique was used for 
constructing the relationship matrix between supplier 
selection criteria and customer requirements, and COPRAS 
method was adopted to prioritize the candidate suppliers. 
Furthermore, a more specific literature review relating to the 
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methods of green supplier estimation and selection can be 
seen in [31, 32]. 


B. APPLICATIONS OF MULTIMOORA METHOD 

In recent years, researchers have made considerable 
extensions of the MULTIMOORA method to handle various 
MCDM problems. For example, Zavadskas et al. [33] put 
forward a modified MULTIMOORA method based on 
IVIFSs to address civil engineering problems. Liu et al. [34] 
extended the MULTIMOORA method under the interval 2- 
tuple linguistic context for evaluating and selecting health- 
care waste treatment technologies. Deliktas and Ustun [35] 
proposed a hybrid method by integrating fuzzy 
MULTIMOORA method and multichoice conic goal 
programming to elect the best student. Hafezalkotob and 
Hafezalkotob [36] developed a MCDM model on the basis of 
MULTIMOORA method and Shannon entropy for managing 
material selection problem. Li [37] developed an extension 
MULTIMOORA method in the context of hesitant fuzzy 
numbers for software selection. Liu et al. [38] applied the 
MULTIMOORA method and fuzzy sets for the prevention of 
infant abduction. Tian et al. [39] developed an improved 
MULTIMOORA' method by integrating simplified 
neutrosophic linguistic normalized weighted Bonferroni 
mean, simplified neutrosophic linguistic normalized 
geometric weighted Bonferroni mean operators and a 
simplified neutrosophic linguistic distance measure. Sahu et 
al. [40] presented an improved MULTIMOORA method in 
interval-valued trapezoidal fuzzy context to evaluate and 
rank computer numerical control machine _ tools. 
Hafezalkotob and Hafezalkotob [41] proposed a target-based 
MULTIMOORA method combined with _ significant 
coefficients for biomaterials selection. Lazauskas et al. [42] 
used AHP, additive ratio assessment (ARAS) and 
MULTIMOORA to assess unfinished residential buildings 
and choose the most appropriate construction project. 
Stanujkic et al. [43] described an extended MULTIMOORA 
by interval-valued triangular fuzzy numbers for selecting 
comminution circuits design scheme. Aytag et al. [44] 
offered a MCDM model based on the MULTIMOORA and 
multi-objective optimization simple ratio analysis 
(MOOSRA) methods to handle the laptop selection problem. 
Gou et al. [45] developed a modified MULTIMOORA 
method based on double hierarchy linguistic term sets to 
select the optimal city in China. 


lil. PRELIMINARIES 
Some basic concepts and operational rules on IVIULSs [13] 
are recalled in this section. 


Definition 1: Let X be a given domain and 5, € S, then an 
IVIULS is denoted by 


A={< x(5,,8; (x),9;(x)) >}, (1) 


where 5, = Ee Sg [38 regarded as an uncertain linguistic 
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variable; 0 (x) and r(x) are the subscripts of the lower limit 


and upper limit to § respectively. The 


x > 
intervals a; :X > D{0, 1] and ¥,:X > D{0,1] respectively 
represent the membership degree and non-membership 
degree of the element x to the uncertain linguistic 
variable S, with the constraint 
0 <sup(i,(x))+sup(,(x))<1,xeX : 

For any element xeX , u,(x) and ¥,(x) are closed 
intervals and their lower points and upper points are 

U 


presented as uz (x) U5 (x) : Vi (x) and ve (x) . Then A can 


be denoted by 
{(((5ueoS0 beh (x),uf (x) [v5 (x).¥4 (x)])) Ire xt, (2) 


where Sy,),5,,, € 5,0 us (x)+ ve (xj, ue (x) =Oand 


L 
vi (x)20. 

For any element x<¢ X , the hesitation interval of the 
element x to the uncertain linguistic 


variable 5, = [Seeaie om is computed as: 


(x) =[ 44 (x)-75 (x) J= [44 (2) 4 (2), 1-05 (2) -5 (2)]- G 


Definition 2° 
Let A= {( (Luss Le (x),u! (x)] [v4 (2).v4 (2)])) re x} 
be an IVIULS, then the 6- 


tuple [54.54 Les (a).u8 (x) DEC) r¥ ()]) is 


called an interval-valued intuitionistic uncertain linguistic 


number (IVIULN). A can be regarded as a collection of 
IVIULNs, thus, 


A={((SqoySey flat (af (x)} Pi ()»4 (9) |xex}- 
Suppose 4, = ( Say Le’ (4).u” (a) Lo" (@), v8 (4 ))) 
and G, = ([ Soca» pL" (4),u” (a,)].[v' (@).v” (@ ))) are 


two IVIULNs and A>0 , the 
of a, and a, are given below [46, 47]: 


1) a, Oa, = ([ Soayeotay>Setarertey >L1-(1-#" (G,))(1-u" (4 )). 


operational rules 


( 
2) 4,4, =([ SyyetaSlayeeta [> ¥* (4, )u* (4) 
[1-(1-v' (4))(1-v' (@,)).1-(1-v" (4) 
3) Aa, = ([ sa) Sita b[ 1 (1a (a) A-(1-w” (a,)) |, 
Lo" (@))’-(0" ad)" J): 
a (Pace iG (4))(u! (a) | 
[1—(1- v4 (4,))'1-(1-»" (4,))']) 
For the purpose of comparing IVIULNs, the expected 
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value and the accuracy function of an IVIULN are 
introduced as follows. 
Definition 3: Suppose 


4=([syay%q bet ade (A) @ 0" (a) isa 


IVIULN, its expected value E (4,) is computed by [13]: 


2 ener (4) 
=a, Liz Ups Liz Urn > 
(9(& )+2(&))x(u" (a) +n” (& )+2—v" (a )-v" (@))/8 
and the accuracy function T (a, ) of a, is defined by 


w(a)eu" (a) v'(a)+o" (a) 
T(a,) = S(a(a)+e(a))/2 | 5) = 2 (5) 


= Sula jeu” (a )+v" (a +0" (&)x(0(a)+2(&))/4° 


Definition 4: Let 
a {[sacay>Ssay |e” (d,),u" (4) ]..v" (4,),v" (a)]), 
a= ([saay Sa) pL (4,),u" (4, )}[v* (4,).v° (4, ))) 


are two IVIULNs. The comparison rules of [VIULNs are 
defined as follows [48]: 
1) IfE(a,)>E(d,), thena, >a,; 
2) If E(a, ) = E(d,), then 
(a) if (4,)>T(G,), thena, >a,; 
(b) ifT (a,)=7 (4, ), then, =4,. 

To aggregate uncertain linguistic information, the 
interval-valued intuitionistic uncertain linguistic weighted 
average (IVIULWA) operator is proposed by Liu [13]. 

Definition 5: Let 
a, = (Be ae (4,),u" (4, )].[v" (4, jj v" (4,)]) 

(i =1, 260) be a collection of IVIULNs, then the 
IVIULWA operator is defined as 


IVIULWA(G,,4,,....4,) = )w,4,, (6) 


j=l 


where w =(W,,W,,.-.5 w,) is the associated weight vector 


of a, (j =1,2,...,2) , which satisfies w, €[0,1] and)" w, =1. 
j=l 


Definition 6: Suppose 
a ([Sotay> Stay} L#” (4,),u” (4,) ].[v" (4).v” (4, )]) ; 
i= ([ Sota Sea Lo” (4,),u" (4) ],[v*(a).v" (a,)]) 


are two IVIULNs. The Hamming distance 
between a, and a, is defined as: 
ies of | (7) 
IV. THE PROPOSED GREEN SUPPLIER SELECTION 
MODEL 


6 +a? (4, )—u" (4,)]+ |r" (@)—v" (4) +r" (@)-»” (@) 
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This section develops a hybrid MCDM approach for 
selecting most appropriate green suppliers under the large 
group environment. In the proposed model, decision makers 
are clustered by employing an ant colony algorithm; 
assessment values of clusters are aggregated by the 
IVIULWA operator; criteria weights are computed by the 
LINMAP method; the alternative suppliers are ranked 
finally based on a modified MULTIMOORA approach. The 
flowchart of the green supplier selection method being 
proposed is shown in Figure 1. 


Phase I 
Cluster decision makers 


Apply ant colony algorithm to cluster decision makers 


== : 


Normalize the cluster decision matrix 


Phase II ' 
Aggregate assessment values of ' i 
clusters 


Utilize the IVIULWA operator to aggregate clusters 


a f 


Provide the united preference relations between 
candidate suppliers 


x 


Define the consistency index and the inconsistency index 


¥ 
Construct a mathematical programming to derive criteria 
weight vector 


- = ; 


Rank alternative suppliers based on the ratio system 


¥ 
’ Rank alternative suppliers based on the reference point 

Phase IV H approach 
Rank the alternative suppliers H 


Phase Ill 
Calculate the weights of criteria 


x 
Rank the alternative suppliers based on the full 
multiplicative form 


Determine the final ranking of the alternative suppliers 


FIGURE 1. Schematic diagram of the proposed model. 
For a green supplier selection problem with m 


alternatives A,(i=1,2,...,m), n criteria OF (JHL2uan)s 
and L decision makers DM, (F122, h> 20) , the 


evaluation matrix by the kth decision maker is denoted 


pk _|[ ~k wk Ea Lou )\; 
as P" = [ Bi 1. , where pj, = (oe Sy Jes Mohn) is 


the IVIULN given by DM, for the alternative supplier A, on 


the criterion C : } with the 
ee L L U U L U 
condition 0<uj;, <1, O0Svy $1, Uiq_ +Vig S1, Ui Si» 
Voi < ve ,5,,5, €S. Next, the procedure of the proposed 
Gijx Gijx 


green supplier selection approach is explained in the 
following subsections. 


B. CLUSTER DECISION MAKERS 

Ant colony algorithm is an effective clustering method 
inspired by the foraging behavior of ant colony. Real ants 
can choose the shortest route between their nest and food 
source. Ants will leave the pheromone in the path and can 
be aware of the existence and strength of pheromone during 
the course of movement. The larger the number of ants pass 
through a given path, the more pheromone the path 
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accumulates [49]. In turn, the number of ants choosing the 
path will increase with a positive feedback effect [50]. 


Finally, the ants will find the shortest path. In this algorithm, 


the clustering center is “food source”. The clustering 
process can be viewed as the process of ants looking for the 
“food source”. Ant colony algorithm is a promising method 
for data classification, which intend to discover a list of 
classification rules (Liang et al., 2016). 

In this paper, we cluster decision makers by using the ant 


colony algorithm. For each supplier A, (i =1, 2,...,m) , the L 
decision makers DM, (k = 1,2,...,L) can be clustered into m 


clustering results based on their evaluations. The algorithm 

under the IVIUL environment is expressed as below: 

1) Transform all the ITVIULNs into expected values as 
input data. 

2) Set parameters: the number of clusters g, the number of 
ants f, the maximum number of iterations f,. , 


pheromone threshold q, and evaporation rate € . 
3) Initialize the pheromone matrix I,,,. The initial value 
is 0.01. 


4) Determine ants’ path according to the value in 
pheromone matrix IT, . If all the pheromone values of 


of each element inI1,,,, 


a sample are less than the pheromone threshold q, then 
the cluster corresponding to the largest pheromone is 
selected. If the pheromone values are greater than the 
pheromone threshold gq, the cluster is determined based 
on the proportion of pheromone in each path of total 
pheromones. 

5) Identify the cluster centers on all criteria of each cluster. 
The cluster center is the average of all samples in the 
cluster with respect to each criterion. 

6) Compute the sum of the Euclidean distances (deviation 
error) of each ant from each sample to its corresponding 
cluster center, and select the path corresponding to the 


smallest deviation error as the best path for this iteration. 


The smaller the deviation error, the better the clustering 
effect. 


7) Update pheromone matrix II,,, .The updated value is 


Lxg 
the original pheromone value multiply (1-«) and plus 
the reciprocal of the minimum deviation error. Choose 
the best path according to the new pheromone matrix 
and perform iteration operations until reach the 
maximum number of iterations f,,,, 


B. AGGREGATE CLUSTERS 
Suppose G* (z =1, 208) is the zth cluster and L, is the 
makers in G* , with the 


number of decision 


g 
condition YL = L,. The evaluation for the supplier A, with 
z=l 


respect to the criterion C; of the cluster G*is the average of 


the assessments given by decision makers, i.e., 
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(8) 


1 1 1 1 
: hex “LE |e mores Dares 


i=1,2,.4,m, J=1,2, 04732 = 1,2, 358 
To eliminate the effects of different dimensions and 
ensure the compatibility of IVIULNs with respect to all 
criteria, we normalize the cluster decision 


matrix P* = [ B; | by using the following formulas: 


ij =([54 E tabla? Pal sili [ «i Uije> uy, | Ly: iz? ve] f eJ,), (9) 
ij = (5c) ae ~ab, fe ea oso Mako) G € J;); (10) 


b= 12554 


where J, and J, denote the sets of benefit criteria and cost 


,m, j =1,2,....n,Z=1,2,...g, 


criteria, respectively, and 


av = max (a/, wot = 12,.. 


jmax 


M,Z = 128) ; 

Generally, the cluster with more decision makers should 
be given a bigger weight. Thus, the weight of each 
cluster G* (z == Per g) can be calculated by 


l 
A.==, z=1,2,...,2. 11 
z g (11) 


Then, we adopt the [VIULWA operator to obtain the 


group normalized decision matrix R = [ aA , where 7, is 


determined by 
i, = IVIULWA(,, Fu?) = Wadi 


& 


Eyes a-fyc-es) ft tTee) 


Si Sa {1 2 
Vig. YAagj. zal kel 
al al 
i=1,2,...5m) J =1,2,...,2) 


where /, is the associated weight of 


£= (are leben bo 12 


the condition A, €[0,1]and A =1. 


C. CALCULATE THE WEIGHTS OF CRITERIA 
The LINMAP is a classical and effective MCDM method 
proposed by Srinivasan and Shocker [51]. It is can be used 
for computing the objective criteria weights by establishing 
a mathematical programming [52]. Therefore, in this study, 
we extend the LINMAP method to the IVIUL context for 
deriving the weight vector of criteria. The specific steps are 
given below: 

Step 1: The decision makers give the united preference 
relations between the candidate 


A,)|A, 2 4,.(2,2=1,2,....m)}. 


The decision makers provide the corresponding pairwise 
comparisons of the candidate suppliers as a whole. The 


suppliers Q = {(A 
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priority relations between alternatives are determined by 
overall judgement rather than on each criterion. 

Step 2: Define the consistency index and inconsistency 
index between objective ranking order and subjective 
preference relation. 

The IVIULN ideal solution (reference point) on each 


re * x * 
criterion are represented as roa (7; Ty yeeeT, ) > 


OO Lou L vu - F 
where r= ess hla ju; ].[¥} Vj Ne = 1, Qe sis, n) is 
the best rating on the criterion C,;. Then the Hamming 


distance between each alternative of (A,,A,) €Q and ideal 


solution r’ is calculated by 


D, = iw, (iF; ) (13) 


p= Swa(s.t) as 


where w = (w,, W 5.25 W, ) is the criteria weight vector. 

For each pair of the alternative suppliers (A,, A, ) € Q, the 
decision makers prefer supplier A, to A, or make no 
difference between A, and A, . Thus, if D, > D,, the ranking 


order determined by D, and D, is consistent with the 


Otherwise, 
determined 


preference relation by decision makers. 
if D.<D, , then the ranking order 
by D. and D, is inconsistent with the preference relation by 


decision makers. Accordingly, the 


indexes (D,-D,) and(D,—D,) are respectively used to 
measure consistency and inconsistency between the 
objective ranking orders determined by D, , D, and 
subjective preference relation [53, 54]. 

(D,-D,) =max{0,(D,-D,)}, (15) 

(D,-D,) =max{0,(D, -D,)}. (16) 
Hence, the total consistency index and inconsistency index 
of all pairs of suppliers can be calculated by the following 
equations: 


H= > (D,-D,) = > max{0,D,-D,}, (17) 


(h,i)eQ (hi)eQ 
B= )° (D,-D,) = > max{0,D, -D}. (18) 
(h,i)eQ (h,i)eQ 
Step 3: Construct a mathematical programming to derive 
the criteria weight vector. 
Let Z,, = max {0, D, — D,} for each pair(A,,A,) €Q, with 
the condition Z,, 20. Then Z,, = D, — D,is obtained. Since 


the total inconsistency index B reflects the group 
inconsistency between the objective ranking order and 
subjective preference relations, then B should be minimized. 
Therefore, we can construct a mathematical programming 
as follows: 


hi 
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Z,, + D, —D, = 0, (19) 
S.t.4Z,, 2 0, 


where 4s > O is used to ensure the total consistency index H 


bigger or equal to the inconsistency index B, and dé >Ocan 
ensure the weight of each crterion greater than 0. 


D. RANK THE ALTERNATIVE SUPPLIERS 
After obtaining the group normalized evaluations of 
suppliers in the second stage, we utilize a modified 
MULTIMOORA method to rank the alternative suppliers in 
this subsection. The ranking process is presented as below. 
Step 1: Rank alternative suppliers based on the ratio 
system approach. 
The evaluation values of ratio system is obtained by 
adding the normalized ratings on all the criteria of each 
supplier, i.e., 


Ee Wi Telia, (20) 


where 7, represents the overall evaluation value of A, with 
regard to all the criteria. The ranking of the alternative 
suppliers are determined by the values Ii =1,2,...,71) in 


descending order. 

Step 2: Rank alternative suppliers by the reference point 
approach. 

In the third phase, we have determined the ieal solution 
as the best rating. Then the weighted distance of each 
supplier to the ieal solution is computed by 


d, =) w,d(F,.F), i=1,2,...,m. (21) 
j=l 


The alternative suppliers are ranked based on 


distances d, (i = 1,2,...,m) in increasing order. 


Step 3: Rank alternative suppliers based on the full 
multiplicative form approach. 

The overall utility of each supplier can be represented as 
an IVIULN by 


U, =| [w,f,, i=1,2,...,m. (22) 


The green suppliers are ranked according to the 


values U, (i =1,2,..., m) in descending order. 
Step 4: Determine the final ranking of alternative 
suppliers. 


By utilizing the theory of dominance, we can integrate 
the three rankings acquired by the ratio system, the 
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reference point and the full multiplicative form to 
determine the final ranking of the m green suppliers. 


V. CASE STUDY 


A. IMPLEMENTATION 
This section applies the proposed model to select the 
optimal timber green supplier for a real estate company in 
Shanghai, China. The related research data show that the 
carbon produced by the real estate industry in China 
accounted for about eight percent of global carbon 
emissions. As facing the great opportunity of green 
tranformation, Alashan SEE ecological association, China 
urban real estate developers strategic alliance, all real estate 
association, Vanke enterprise and Landsea green real estate 
co-sponsored the action of green supply chain in the real 
estate industry on June 5, 2016. The white lists of qualtified 
suppliers are determined by real estate companies 
participating in the green supply chain action according to 
green standards, which include steel suppliers, cement 
suppliers, timber suppliers and aluminum suppliers. This 
study aims to assisst the real estate company to find out the 
most appropriate timber supplier. 

After initial screening, four suppliers (A,,i =1, 2,...54) are 


remained for further assessment and selection. Based on a 
literature review and expert interviews, five evaluation 
criteria (C, J=h Bia) , product quality, purchase cost, 
technology capability, green degree and delivery level 


(Liao, Fu, & Wu, 2016), are identified for the green 
supplier selection. Besides, an expert group consisting of 


twenty-two decision makers, (DM,,k =1,2,...,22) , is 
invited to evaluation the performance of the four alternative 


suppliers with respect to each criterion. The linguistic terms 
set S is used in the performance evaluation of the suppliers. 


S= {sy =very poor, S, = poor, Ss, = medium poor, s, = fair, 
S,=medium good, s, = good, Ss, =very good}. 


The assessments of the four alternative suppliers on the 
five crieria provided by all the decision makers are shown 
in Tables 1-4. First, we transform the IVIULNs into their 
expected values as input data. Then the parameters of the 
ant colony algorithm are set as: the number of 
clusters g =5, the number of ants f = 1000, the maximum 


trax = L000 , 
threshold g =0.9 , and evaporation rate ¢ =0.1. By using 


number of _ iterations pheromone 


the ant colony algorithm, the clustering results of decision 
makers are shown in Table 5. 

Based on Equation (8), the evaluation values for the four 
suppliers of the five clusters are calculated as shown in 
Table 6. Next the evaluation values of the five clusters are 
normalize by utilizing Equations (9)-(10), and the results 
are shown in Table 7. By Equation (12), the normalized 
evaluations of the five clusters are aggregated into the 


collective normalized evaluation matrix R=|F, |, 5 as 
x, 


shown in Table 8. 
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The united preference relations of the candidate suppliers 
given by the twenty-two decision makers 


are: Q= {(1,3),(2.1),(2, 4),(1,4)}. The ideal solution is set 
as ([s,, CAs 1],[0, 0}) for each criterion. Then the following 


linear programming is eatablished to calculate the weights 
of criteria: 

min Z,,+Z,,+Z,,+Zy4 

0.839w, —0.768w, + 2.792w,+1.117w, +0.924w, 21, 

0.019w, +0.229w, +0.574w, +0.057w, +0.130Ww,+Z,, > 0, 

0.252w, —0.664w, + 0.537, + 0.231w, + 0.133, + Z,, 2 0, 

St.) 0.273w, —0.282w, +0.553w, + 0.386w, + 0.265w, + Z,, 2 0, 

0.294w, —0.052w, +1.127w, + 0.443w, + 0.395w, + Z,, = 0, 

w, 2 0.01,(j =1,2,...,5), 


w,+w,+w;+w,+w, =1. 


By solving the above model, the weights of the five criteria 
is derived as w = (0.28,0.10,0.10,0.42,0.10) . 

In the fourth stage, the modified MULTIMOORA 
method is used to rank the four timber green suppliers. 
Firstly, based on the ratio system approach, the ranking 
order of the alternative suplliers are determined by 
Equation (20), and the result is: 


7, = ([Soes7> Sos0s ]-[0-608, 0.742], [0.099, 0.206] 
7 = ([5o098>Sosio |-[0-910, 0.953] ,[0.000, 0.000] 
Fs = ([ S507» Sore |-[0.599, 0.734] ,[0.122, 0.229] 
T= ([sosss+So2s | [0.576, 0.717], [0.146, 0.241). 


Based on the comparison rules of [VIULNs in Definition 4, 
the ranking of the four suppliers is derived as A2> A,> A3> 
A 


\, 
F ): 
; : 


Secondly, in accordance with the reference point 
approach, the weighted distance of each supplier to the 
ideal solution is computed by Equation (21), and the result 
is: 
d, =1.581,d, =1.458,d, =1.750,d, =1.874. 

Thus, the four timber supplilers are ranked as A2> A1> A3> 
Ag. 

Thirdly, according to the the full multiplicative form 
approach, the overall utility value of each supplier is 
determined by Equation (22). 


D, = ([Scoesn0 Srapnes [8.472*10, 6.566*10~ ],[0.999,0.999 


se 


S S : 
[ 1.168*107 ? "2.94710 


i }): 
( ].[9.033*10°,4.075*10~ ],[0.997,0.999]), 
Ts = ([ 5, e065, apreos [1-834 *10°,9.952*10° ],[0.999,0.999]}, 
( 5s nueio® |>[5-696*10 *, 5.342*10* ],[0.999,0.999]). 


[Sewers 2"'1.716*10~ 
As a reslult, the ranking of the four suppliers is derived as 
A2> A\> A3> Aa. 

To summarize, the final ranking of the four timber 
suppliers is A2> A;> A3> A4 based on the dominance theory. 
Therefore, Az is the optimum green supplier among the four 
alternatives for the considered application. 
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TABLE | 
PERFORMANCE VALUES OF THE FIRST SUPPLIER Aj BY THE DECISION MAKERS 
Decision _ Criteria 
makers C C2 C3 Cy Cs 
DM: ([ss.,].[0.7,0.8],[0.1,0.2]) ([s,.85],[0.6,0.7],[0.1,0.2]) _([;,5,].[0.8,0.9].[0.1,0.1]) ([s,.5, ],[0-5,0.6],[0.1,0.3]) ([ss., ].[0-8,0.8],[0.1,0.2]) 
DM? ([s,.s, ].[0.4,0.5],[0.2,0.3]) ([s,.85],[0-6,0.8],[0.1,0.2]) —_([s,,, ],[0.5,0.7], [0.1,0.2]) ([s,.85],[0-6,0.8],[0.0,0.2]) _([s,,,],[0.7,0.7].[0.1,0.2]) 
DM; ([s,.85].[0.6,0.6],[0.1,0.3]) ([s5.56],[0.7,0.8],[0.1,0.2]) —_([s,,8,].[0.7,0.9] ,[0.0,0.1]) ([s,.8,],[0.6,0.7],[0.0,0.2]) —_([s,,,],[0.5,0.6],[0.1,0.2]) 
DMn die 5,],[0.7,0.9], [0.1,0.1]) ([sussh0.7,09}0. 1,0.1]) ((sss¢}{0.7,038,[0.1,0.1) ([ss}{0.6.0.7}[0.1,0.3) ((sees,}0.6,0.7] 0.1.0.3) 
DM» ([s,-8.].[0.7.0.8],[0.1,0.2]) ([5s.5.],[0-7,0.8],[0.0,0.1]) —_([ 5,5, ],[0.6,0.7], [0.1,0.2]) ([54.5,],[0-7,0.7],[0.1,0.2]) _ ([s,,,],[0.7,0.9],[0.1,0.1]) 
TABLE 2 
PERFORMANCE VALUES OF THE SECOND SUPPLIER A? BY THE DECISION MAKERS 
Decision _ Criteria 
makers Cc; C2 C3 Cy Cs 
DM: ([s..5].[0-6,0.8],[0.1,0.2]) ([s,.55].[0-5,0.7],[0.1,0.3]) _([s,,,].[0.8,0.8],[0.1,0.2]) ([s,.8, ],[0.5,0.6],[0.1,0.2]) ([s,.55].[0-6,0.8],[0.1,0.2]) 
DM? ([s,.5].[0.7,0.8],[0.2,0.2]) ([s,.85],[0.6,0.7],[0.1,0.2]) (Ls... ].[0.7,0.7], [0.1,0.2]) ([s;.s, ],[0.5,0.7],[0.0,0.2]) ([s,.55].[0-5,0.6],[0.2,0.3]) 
DM; ([ss.8, ].[0.6,0.8],[0.1,0.2]) ([s,.55].[0.6,0.7],[0.2,0.3]) _([,,,].[0.8,0.9],[0.0,0.1]) ([s,.8,].[0-5,0.7],[0.1,0.2]) ([s,.5s ].[0-8,0.9],[0.1,0.1]) 
DMn, (sws.){0.7,038},[0.1,02]) eal [0.6,0.9],[0.0,0.1]) ( s,,8.],[0.6,0.7],[0.1,0.1]) (ss s,],[0.8,0.8].[0.1,0.2]) (se s,],[0.7,0.7],[0.1,0.3]) 
DM» ([s,.55].[0.7,0.8],[0.1,0.1]) ([s,.5,].[0-8,0.9],[0.0,0.1]) ([5..5,] [0-6,0.8],[0.1,0.2]) ([s,.8, ],[0.5,0.7],[0.1,0.3]) ([s..5.].[0-6,0.8],[0.0,0.2]) 
TABLE 3 
PERFORMANCE VALUES OF THE THIRD SUPPLIERA3 BY THE DECISION MAKERS 
Decision Criteria 
makers Ci C2 C3 Cy Cs 
DM: ([s,.55].[0.6,0.7],[0.1,0.3]) {5,5 ],[0.5,0.6],[0.2,0.3]) —([.s,,8,].[0.7,0.8],[0.1,0.2]) —_([,,55].[0.7,0.8],[0.1,0.2]) _([,,, ].[0.6,0.7],[0.1,0.2]) 
DM2 ([5.5.],[0-7,0.8],[0.1,0.2])  ([s,,5,],[0-8,0.9],[0.1,0.1]) — ([s5,5].[0.6,0.7],[0.1,0.2]) _([,,,].[0.5,0.7],[0.1,0.2]) (s,s, ],[0.5,0.6],[0.3,0.4]) 
DM3 ([s,.55].[0.7,0.8],[0.1,0.2]) — ([s,.s, ].[0.6,0.7],[0.2,0.2]) —_({s,,s,].[0.7,0.8].[0.0,0.1]) _([s,,,].[0.6,0.8],[0.1.0.2]) —_([s,,,],[0.8,0.9],[0.0,0.1]) 
DMay ((sus,)[0.5:0.6},[0.1,0.2) ( 8,,8,],[0.6,0.8], [0.0,0.1]) a s;,5,],[0.6,0.7],[0.2,0.2]) ((sses,}0.7.08},[0 1,0.2}) (ses 06,0.7] 0.1.0.3) 
DM ([s,.55].[0-8,0.9],[0.1,0.1]) ([ss.85].[0-8,0.9],[0.0,0.1])  ([s,,8,],[0.7,0.8],[0.1,0.2]) — ([,,s,].[0.5,0.7],[0.1,0.2]) _([s,,s,],[0.6,0.8],[0.1,0.2]) 
TABLE 4 
PERFORMANCE VALUES OF THE FOURTH SUPPLIERA4 BY THE DECISION MAKERS 
Decision Criteria 
makers Ci Cr C3 Cy Cs 
DM: ([s5.56].[0-6,0.6],[0.1,0.3]) — ({s,,8,].[0.5,0.6],[0.2,0.4]) —_([s,,s,],[0.7,0.7],[0.1,0.2]) ({s,,8,],[0.6,0.8],[0.1,0.2]) —_([,,, ].[0.6,0.7], [0.1,0.3]) 
DM, ([s,.55].[0-8,0.8],[0.1,0.2]) — ({s,,8,],[0.6,0.7],[0.0,0.1]) —_([s,,,],[0.5,0.7],[0.2,0.2]) ([s,,8,],[0.5,0.6],[0.1,0.3]) —__([s,,, ].[0.6,0.7],[0.3,0.3]) 
DMs3 ([s,.8,].[0.7,0.9],[0.1.0.1])  ({s,,s,].[0.6,0.8],[0.1,0.2]) _([s,,,],[0.7,0.8],[0.0,0.2]) ——([s,,8,],[0.5,0.6],[0.1,0.3]) __[s,,8,].[0.7,0.7], [0.0,0.1]) 
DM (sus. )0.5.0.6]},[0.3,0.3]) ([su.8}0.6.08},[0.1,02]) ((sssh{06.0.7],0.2,03}) (sss) {0.7,0.7],[0.1,02)) ((ss,}{06,09},[0.1,0.1) 
DM2. 


([s,.85].[0.6,0.8],[0.0,0.1]) 


([ss.s].[0.6,0.8],[0.0,0.2]) 


([s,.55].[0.7,0.7], [0.1,0.2]) 


([s;.8, ].[0-5,0.7],[0.2,0.3]) 


([s,.5,].[0.5,0.7].[0.2,0.3]) 
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TABLE 5 


CLUSTERING RESULTS OF THE FOUR SUPPLIERS. 


Suppliers Clusters 


G! CG G G* G 

A DM,, DM,, DM, DM 0, DM;, DMs3, DM,, DM, DM3, DM 
DM;, DMo, DM 4, DMs, DM,3, DM, DM, 

Ad DM,, DM3, DM,, DM;, DM, DM, DM, DM,3, DM» 
DM;, DMs, DMs3, DMo, DM\4, DM, DM 
DM6, DMs, DM DMs 

A; DM,, DM3, DM,, DMs, DM,, DMs, DM;3, DM\a, DM, 
DMg, DM,2, DMo, DM 0, DM 3, DM, DM», DM, 

A4 DM,, DM3, DM,, DM;, DM;3, DM, DMs, DM, DMs, DMs, 
DMs;, DMg, DMo, DMo, DM\4, DM DM), DM 
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TABLE 6 


EVALUATION VALUES BY THE FIVE CLUSTERS 


Suppliers Clusters _ Criteria 
Ci Cy C3 Cy Cc. 5 


Aj G' ([ss.s2+8ssn ]-[0.671, 0.829], [0.086,0.157]) ([53457-54n4 ].[0-586, 0.729], [0.114,0.186]) ([sass7 +9557 ]-[0.714, 0.842], [0.071,0.129]) [5.000 +424 ],[0-600, 0.757], [0.071,0.200]) ([51000 +4214 ]-[0-629, 0.757], [0.114,0.200}) 
C ([5s500-84233 }-[0-550, 0.700] ,[0.150,0.233]) ([53s33+ 84332 ]-[0.633, 0.767], [0.083,0.183]) ([5ss33> S483 ] -[0.567, 0.783], [0.133,0.183]) ([55.500+S4167 ]-[0-583, 0.733] ,[0.050,0.200]) ([Sa67 +5167 ]-[0-650,0.783],[0.100,0.167]) 
G [5500095250 }[0-650, 0.800] ,[0.100,0.150]) ([5,250» 5s sw ]-[0-650, 0.800], [0.075,0.175]) [5.504250 ]-[0-600, 0.800], [0.075,0.175]) [55250 %4 00 ]-[0-675, 0.775] ,[0.100,0.150]) ([55250 4500 }-[0.575, 0.725], [0.125, 0.225]) 
G* [5,000 S333 ]-[0-567, 0.767] ,[0.100,0.233]) [5.000 8607 ]-[0-600, 0.733], [0.100,0.200]) [5,607 5533s ]-[0-767, 0.867] ,[0.067,0.133]) ([5.33>%000 } [0.667, 0.767] ,[0.067,0.167]) ([4000» %4a00 ]+[0-567, 0.700], [0.100,0.200}) 
G [5,000 sco ]-[0-650, 0.700], [0.100,0.250]) [55.000 Sscco }-[0-700, 0.800], [0.050,0.150]) ([55000 55500 }-[0-650, 0.800]_,[0.050,0.150]) ([5s500> ao ]-[0-650, 0.700], [0.050,0.200]) (E5000 Ss. ]>[9-600, 0.750], [0.100,0.150]) 
Az G' ([1250 550 ]-[0-663, 0.788], [0.125,0.175]) ([53s75+ 84500 }-[0.575, 0.713], [0.100,0.213]) ([a:25+5000 }-[0-688, 0.800], [0.088,0.150]) ([53.500- 84.125 } [0-575, 0.725], [0.088, 0.200]) ([.000 47s ]-[0-663, 0.750], [0.125,0.225]) 
CG ([5000 500 }-[0.580, 0.780], [0.120,0.200}) ([5s.200Sscco ]-[0-700, 0.840], [0.080,0.140]) (L500 85.400 }-[0.640, 0.780], [0.100,0.180]) ([5.s00>%s 00 ]+[0-600, 0.780]. [0.080,0.200]) ([.000 +5000 ][0-650, 0.700] ,[0.100, 0.250}) 
G ([5ss00 +%4 sco ]-[0-600, 0.760]. [0.120,0.220]) ([3.00 5300 } [0-620, 0.760], [0.120, 0.200}) ([5.200 sao ]-[0-680, 0.780], [0.140,0.180]) ([52200 2:00 ][0-660, 0.760] ,[0.080,0.180}) ([2s00 +300 ]-[0-660, 0.780] ,[0.100, 0.200}) 
G [5465755333 }-[0-733, 0.800] ,[0.067,0.200]) [5.461 5s.07 ][0-667, 0.800], [0.067,0.167]) [55000 aco ]-[0-600, 0.733] ,[0.100,0.200]) [5,000 8607 ][0-667, 0.833], [0.100,0.100}) (L5s000> 83467 ]+[0-600, 0.733], [0.100,0.267]) 
G [55.000 %saco ]-[0-500, 0.700], [0.100,0.100]) [55000 4a ]+[0-600, 0.900], [0.100,0.100]) [5,000 %5cao ]-[0-800, 0.900], [0.000,0.100]) (E5000 %s.u0 }-[0-800, 0.900] ,[0.100,0.100]) (E5500 Sea ]+[0-800, 0.800] ,[0.100,0.200}) 
A3 G' ([5.167°5so0 ]+[0-617, 0.700], [0.150,0.267]) [55667 84,167 }[0-633, 0.717], [0.117,0.200}) [54533595333 ] [0.650, 0.767]_,[0.083, 0.183]) (Ls.167> 5.0 ]-[0-633, 0.800] ,[0.083,0.183]) ([ss1s1 84167 ]-[0-600, 0.667], [0.133, 0.283]) 
CG ([s.000 +667 ]-[0-667, 0.783], [0.083, 0.150]) ([3667 4667 | [0-550, 0.717], [0.117, 0.267]) ([532:3+ 8382 ]-[0-600, 0.750], [0.100, 0.200]) ([53323+ 84000 } [0-600, 0.733] ,[0.167,0.267]) ([51333- 8533 ]-[0.650, 0.800],[0.067,0.167]) 
G ([53250 1250 ]-[0-675, 0.750], [0.100,0.175]) ([53250 +4500 ][0-625, 0.800], [0.075,0.150]) ([52250 +3250 ]-[0-625, 0.750], [0.100,0.200}) ([52250 53250 ]-[0-650, 0.775], [0.100,0.200}) ([52250+53250 ][0-625, 0.775], [0.075, 0.200]) 
G ([s.200 +5200 }-[0-660,0.760],[0.120,0.220]) (L500 500 } [0.700, 0.800] ,[0.060,0.140}) ([s.200 +5200 ]-[0-640, 0.760], [0.100,0.220]) ([s.400 5.400 }-[0-600,0.760],[0.100,0.180]) (L000 +400 ]-[0-600, 0.760] ,[0.120, 0.240]) 
G [5.000 Sa ]-[0-500, 0.600], [0.200,0.300]) [55000 8so00 }-[0-700, 0.900],[0.000,0.100]) [5.000 82.9 ][9-600, 0.700], [0.100,0.300}) [5.000 Ssc00 ]-[0-500, 0.800], [0.100,0.200]) ([5.cco +2000 | >[0-700, 0.700], [0.100,0.300]) 
A4 G' ([5,3332 9523s ]-[0.567, 0.700] ,[0.133, 0.267]) [5.167 54s00 }[0-533, 0.633] ,[0.100, 0.267]) [5.167 S167 }-[0-583, 0.700] [0.133,0.217]) ([s.107 5483s ] -[0-617, 0.767], [0.083, 0.200}) ([52533+ 5:33 ]-[0.633, 0.700], [0.150,0.267]) 
C ([54533»%so00 }-[0.650, 0.800]  [0.133,0.183]) ([5.433+ 9:33 ]-[0.583, 0.717], [0.150,0.217]) [5,500 95:35 }-[0-650, 0.767], [0.150,0.200]) ([55500»%4 eo ]-[0-667, 0.767], [0.117,0.217]) ([54.333+so0 ] >[0.583, 0.733], [0.100,0.217]) 
G ([3250 +1250 }-[0-600, 0.825], [0.100,0.150}) (L000 5.000 ]-[0-600, 0.825],[0.125,0.175]) ([2250 53250 ].[0-725, 0.800], [0.050,0.175]) ([5:250+8.250 ]-[0-600, 0.725], [0.125,0.225]) ([2.000 3.000 ][0-625, 0.800], [0.050,0.125]) 
G ([2467 53667 ]-[0-500, 0.633] ,[0.233,0.333]) [52233-53229 ]-[0.633, 0.833], [0.033,0.133]) ([3000 83.607 }.[0.567, 0.667] ,[0.167,0.267]) ([5s000 53233 ]-[0.600,0.667],[0.167,0.200]) ([5.533» 332s ]-[0-533, 0.800], [0.133, 0.200]) 
G ([8sa33+8sax9 ]-[0-633, 0.767],[0.100,0.167]) ([.067 +5233 }-[0-600, 0.733], [0.033, 0.200]) ([3667 51.67 ] [0-633, 0.800], [0.133,0.167]) ([5s000+%saao ]-[0-533, 0.667], [0.200,0.300]) [100054333 ]-[0-567, 0.767], [0.167,0.233]) 
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TABLE 7 


NORMALIZED EVALUATION VALUES BY THE FIVE CLUSTERS 


Suppliers Clusters _ Criteria 
Ci QC C3 Cy Cs 


Aj G' ([soss7 061 ]-[0.671, 0.829], [0.086,0.157]) ([Soars-S0as7 ].[0-114, 0.186], [0.586,0.729]) ([50s72- 51.000 ]-[0.714, 0.843], [0.071,0.129]) ([50467 S076 ]-[0-600, 0.757] ,[0.071,0.200]) ([o.667 Sos ]-[0-629, 0.757], [0.114,0.200]) 
C ([Sosos> Sour ].[0-550, 0.700], [0.150,0.233]) ([Soars S026: ]-[0.083, 0.183], [0.633,0.767]) ([Sosss-Sossr ]-[0.567,0.783] ,[0.133,0.183]) ([so.ss3> 060s ]-[0.583, 0.733], [0.050, 0.200]) ([o.04 Sosa: ]-[0-650, 0.783], [0.100,0.167]) 
G ([5.502 S090: | [0.650, 0.800], [0.100,0.150]) [5008350208 }-[0.075, 0.175], [0.650, 0.800]) ([50.625>0as2 }-[0.600, 0.750], [0.125,0.250]) [50625 89730 }-[0-675, 0.775] ,[0.100,0.150]) ([0625 0750 ]»[0.575, 0.725], [0.125,0.225]) 
G* [5.680 90247 ]-[0-567, 0.767], [0.100,0.233]) ([50222» Sass }-[0-100, 0.200], [0.600, 0.733]) ([50s36 S087 ]-[0-767, 0.867] ,[0.067,0.133]) ([soss0 51000 | [0-667,0.767],[0.067,0.167}) ([0.667 Sae67 }+[0-567, 0.700], [0.100,0.200]) 
G [5.680 Sose2 |+[0-650, 0.700] [0.100,0.250]) ([So.00 $0167 }-[0-050, 0.150], [0.700, 0.800]) ([So07>Soos7 ]+[0-650, 0.800], [0.050,0.150}) ([50.553» 0667 }-[0-650, 0.700]. [0.050,0.200]) ([50500> S607 }[0-600, 0.750], [0.100,0.150]) 
Az G' ([sos10+S0ous }-[0-663, 0.788], [0.125,0.175]) ([S0as0-S0ass ]-[0-100, 0.212], [0.575,0.713]) ([S0240 +007 ]-[0-688, 0.800] ,[0.088, 0.150]) ([Sossa> Sess ]-[0.575, 0.725], [0.088, 0.200]) ([So067>Sos1s ]-[0-663, 0.750], [0.125,0.225]) 
CG ([Sos62 51.000 ]-[0.580, 0.780], [0.120,0.200]) ([50000> 0.1333 }-[0-080, 0.140],[0.700,0.840]) ([0s07 Soe ]-[0-640, 0.780], [0.100,0.180]) ([50633 067 ]-[0-560, 0.700]. [0.140,0.280]) ([Sasu0  Sas67 ]+[0-600, 0.780] ,[0.080, 0.200]) 
G ([o4ss Sosa }-[0-600, 0.760], [0.120,0.220]) ([5o.67>Soass ]-[0-120, 0.200], [0.620, 0.760]) ([50.s24>S0218 }-[0-680, 0.780], [0.140,0.180]) ([S067 S067 ][0-660, 0.760] ,[0.080,0.180}) ([ 0.467 +0433 ]-[0-660, 0.780], [0.100,0.200]) 
G [5040 S019 }-[0-733,0.800] ,[0.067,0.200]) ([50.50 50.56 }-[0-067, 0.167], [0.667, 0.800]) ([50.538» 018 ]-[0.600, 0.733], [0.100, 0.200]) [5.607 5077s ]-[0-667,0.833],[0.100,0.100]) (L500 S611} [0.600, 0.733]. [0.100, 0.267]) 
G ([50s12 Seo }+[0-500, 0.700] ,[0.100,0.100}) ([50.33> 5000 }-[0-100, 0.100],[0.600, 0.900]) [50.218 S0s07 ]-[0-800, 0.900], [0.000,0.100]) [504335033 ]-[0.800, 0.900], [0.100,0.100]) ([50533> 51.000 |»[0-800, 0.800], [0.100, 0.200]) 
A3 G' ([5oris>Sose2 }-[0-617,0.700]  [0.150,0.267]) ([50.06> Sosa }-[0-117,0.200]_ [0.633,0.717]) ([50.798» 0957 }-[0-650, 0.767], [0.083, 0.183]) ([5.604> Sosss }-[0-633, 0.800], [0.083, 0.183]) ([5o.so8» Soe0s |+[0-600, 0.667], [0.133,0.283]) 
CG ([so4s0- Sosos } -[0-667, 0.783] ,[0.083,0.150]) ([50222+S0ass }-[0-117, 0.267], [0.550,0.717]) ([so.sox -Soass } [0-600, 0.750], [0.100,0.200]) ([50.ss6- S007 ]-[0-600, 0.733] ,[0.167,0.267]) ([5o222+Sosso ]-[0-650, 0.800] ,[0.067,0.167]) 
G ([So17 +0732 ].[0-675, 0.750], [0.100,0.175]) [50250-0375 ].[0-075, 0.150], [0.625,0.800]) ([s.403-S0ss3 ]-[0.625, 0.750], [0.100, 0.200]) ([Soss 002s ] [0-650,0.775],[0.100,0.200]) ([Sosrs+Sosu2]-[0-625, 0.775], [0.075, 0.200]) 
G ([Sor21-Sos00 ]-[0-660,0.760],[0.120,0.220]) ([50033 S000 } [0-060, 0.140], [0.700, 0.800}) ([50s1-S0oxs } [0-640, 0.760], [0.100,0.220]) ([50233 50000 }[0-600, 0.760], [0.100,0.180}) [50.67 Saco ];[0-600, 0.760] ,[0.120,0.240}) 
G [50345 S017 ]-[0-500, 0.600], [0.200,0.300]) ([50.500» 50 }+[0-000, 0.100], [0.700,0.900]) ([50.350» 99339 ]-[0-600, 0.700], [0.100,0.300]) [505335000 }-[0.500,0.800]_,[0.100,0.200]) ([0167 S033 }-[0.700, 0.700], [0.100, 0.300]) 
A4 G' ([son41 Sooi9 ]-[0-567, 0.700], [0.133,0.267]) [502502 50306 |-[0-100, 0.267], [0.533,0.633]) ([50.s6s» 50248 ] -[0.583, 0.700], [0.133,0.217]) ([5.604> Sosa |+[0-617, 0.767] ,[0.083,0.200]) ([50.472» 0630 }-[0-633, 0.700], [0.150,0.267]) 
C [50.47 50302 |+[0-650, 0.800], [0.133,0.183]) ([5.361+ 050s ]>[0.150, 0.217], [0.583,0.717]) ([5os0s » 50957 ]-[0-650, 0.767], [0.150,0.200}) ([5o.ss3» 0750 }-[0-667,0.767],[0.117,0.217]) ([50222>Sasss ]-[0-583, 0.733], [0.100,0.217]) 
G ([0.s00 0233 } [0-600, 0.825],[0.100,0.150]) ([So167 +0333 }-[0-125, 0.175], [0.600, 0.825]) ([80.104 +0473 ]-[0-725, 0.800], [0.050,0.175]) ([0.s12 5000s }-[0-600,0.725],[0.125,0.225]) ([5o233-5.500 ][0-625,0.800]_,[0.050,0.125]) 
G ([80.s0 +0432 }-[0.500, 0.633], [0.233, 0.333]) ([S0.04 S011 ]-[0.033, 0.133], [0.633,0.833]) ([50.s38-So4ss }-[0.567, 0.667], [0.167,0.267]) ([50.s00 $0556 ]-[0-600, 0.667], [0.167,0.200]) ([Soss0- 0556 ]-[0-533, 0.800], [0.133,0.200]) 
G ([So47+Soo19 }-[0-633,0.767]_,[0.100,0.167]) ([50.11+ 50.222 ]-[0.033, 0.200], [0.600, 0.733]) ([So4ss>Sosas ] -[0-633, 0.800], [0.133,0.167]) ([0.s00> Soe67 ]+[0-533, 0.667], [0.200, 0.300]) [80057 > 50222 ],[0-567, 0.767], [0.167,0.233]) 
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TABLE 8 


THE NORMALIZED COLLECTIVE EVALUATION MATRIX 


Suppliers C; C2 

Ai ([S0.071»00s | +[ 0-639, 0.783]. [0.098, 0.185]) ([ 50.196» So3ie |-[9.092, 0.183] ,[0.619, 0.756]) 
As ([ 50.205» 50932] -[0.639, 0.781], [0.105,0.191]) ([ 50524» S430 }-[0.099, 0.179] ,[0.621,0.779]) 
A3 ([ 50.682» Sasi7 }-[0.649, 0.747], [0.113,0.201]) ([0255»0350 ]»[0.222, 0.330] ,[0.446, 0.637]) 
As ([o.s2> S030 ]>[0.606, 0.761] ,[0.128, 0.202]) ([ 003» Suz ]>[0.091, 0.196], [0.591,0.741]) 
Suppliers C3 Cy 

Ai ([Saso4>Sooss }-[0-694, 0.829] ,[0.075,0.145]) ([ 50115024 ]-[0.639, 0.756] ,[0.068, 0.179]) 
Ap ([ 50.261» as01 ]-[1-000, 1.000], [0.000, 0.000) ([50102»0779 ]-[0.693, 0.826], [0.112,0.149]) 
As ([.50.650»079s ] [0-630 0.756], [0.095,0.204]) ([ 0.614» 0761 ]»[0-615, 0.771], [0.107,0.208]) 
Ag ([So614>Sosor ]+[0-635, 0.759] ,[0.124,0.197]) ([ 50583» 5020s ]-[0-616, 0.735], [0.122,0.221]) 
Suppliers Cs 

Al ([ 0653+ 0776] >[0.616, 0.753] ,[0.108, 0.187]) 

As ([025»oss2] [0.707, 0.779], [0.174, 0.285]) 

As ([5o.525» S023 ]+[0.624, 0.754] ,[0.095, 0.221]) 

Ag ([0511> 50402] [0-576, 0.758] ,[0.134,0.225]) 


B. COMPARISONS AND DISCUSSION 

To reveal the effectiveness of the developed hybrid 
MCDM model for green supplier selection, a comparative 
study is carried out in this subsection. Based on the same 
case example, the following green supplier selection 
approaches are selected for the comparison analysis: the 
fuzzy TOPSIS [55], the intuitionistic fuzzy GRA (IF-GRA) 
[56], and the interval 2-tuple linguistic VIKOR (ITL- 
VIKOR) [24]. Table 9 displays the results of the four timber 
suppliers according to the listed methods. 


TABLE 9 
RANKING RESULTS OF DIFFERENT METHODS 
Suppliers Fuzzy IF-GRA ITL-VIKOR Proposed 
TOPSIS Method 
Ay 2 2 2 2 
Ay 1 1 1 1 
A3 3 4 4 3 
Ag 4 3 3 4 


From Table 9, it can be observed that the top two green 
suppliers derived by the three comparative models and the 
presented method are exactly the same. Therefore, it 
validates the green supplier selection framework designed 
in this study. Besides, the ranking derived by the proposed 
approach is consist with the one by the fuzzy TOPSIS 
method. The fuzzy TOPSIS method adopted trapezoidal 
fuzzy numbers to evaluate the performance of suppliers. 
However, it is difficult to give approximate ratings to match 
the linguistic phrases in the primary expression domain. 
Moreover, the relative importance of the distances between 
each alternative to the positive ideal solution and the 
negative ideal solution is not considered in the TOPSIS 
method. 

In addition, there are some differences for the rankings 
obtained by the IFS-GRA, the ITL-VIKOR and _ the 
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proposed method. According to the proposed method, A3 
ranks before Aq while the IF-GRA and the ITL-VIKOR 
methods give opposite ranking. The reasons for the 
inconsistence may lie in the weakness of the two methods. 
For the IF-GRA method, IFSs are denoted as crisp values. 
In practice, it is difficult to determine the membership and 
non-membership degrees precisely. In addition, the GRA 
method does not take the grey relation coefficient of 
alternatives to the negative ideal solution into account. For 
the ITL-VIKOR method, the VIKOR method does not 
consider the relative importance of the strategy of 
maximum group utility and the individual regret. 

Based on the above comparison analysis, it can be seen 
that the ranking result derived by the proposed model is 
more precise and credible. Comparing with the listed 
methods, the green supplier selection model proposed in 
this research has the distinctive advantages as follows: 

1) A large group of decision makers are involved to cope 
with the green supplier selection problem, which is 
more practical in the real-life situation. By utilizing the 
ant colony algorithm, large number of decision makers 
can be clustered into several subgroups rapidly and 
effectively. 

2) By applying IVIULSs, the uncertainty and vagueness of 
decision makers’ evaluations can be captured in the 
proposed model. This makes decision makers to express 
their opinions more flexibly and the assessment 
procedure easier to be performed. 

3) The proposed method extends the LINMAP method for 
determining the objective weights of evaluation criteria. 
It is useful to deal with conflicting preference relations 
between alternatives given by decision makers. 

4) The modified MULTIMOORA method under IVIUL 

environment determines the ranking results based on 
three perspectives (i.e., the ratio system, the reference 
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point method and the full multiplicative form). Thus, a 
more precise and reliable ranking of green suppliers 
can be acquired according to this method. 


VI. CONCLUSIONS 

In this study, we developed a hybrid MCDM model to 
address the green supplier selection problem involving 
large group of decision makers. As a new representation of 
uncertain linguistic information, IVIULSs were utilized for 
dealing with decision makers’ diversity assessments of 
suppliers. Ant colony algorithm is utilized to cluster the 
large number of decision makers into subgroups. An 
extension of the classical LINMAP method is used for 
computing the optimal weights of evaluation criteria. 
Finally, a modified MULTIMOORA approach is applied to 
rank the alternative green suppliers. A real-life example of a 
real estate company is implemented to reveal the efficiency 
of the proposed large group green supplier selection 
approach. The results show that the proposed approach can 
better reflect the hesitation and fuzziness of experts’ 
evaluations and obtain a more accurate ranking order of the 
candidate suppliers. 

In this work, the evaluation criteria are supposed to be 
independent in the green supplier selection. Future research 
should be conducted to address the inter-dependence among 
evaluation criteria in GSCM. Besides, it is prospective to 
apply heterogeneous information for evaluating green 
suppliers since suppliers’ performance on different criteria 
may be expressed in different forms. Finally, the developed 
model can be used to handle other large group decision 
making problems, selection, 
resource management and factory location. 


including robot human 
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